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^ \ High precision measurements of the fusion excitation functions for the reactions ^°Ca 

+ ^^^Pt, ^^^Os clearly demonstrate that projectile excitation significantly modifies the 
^ . potential barrier distribution. In sharp contrast, fusion of ^^O + ^^^Sm appears to show 

c ■ 

no influence of the projectile excitation on the shape of the barrier distribution. These 

><' 
Vh ' apparently conflicting conclusions are reconciled in this work, using realistic coupled- 

c^ \ 

channels calculations, which show that high energy states produce an adiabatic potential 

renormalisation. This result indicates that adiabatic effects restrict, in a natural way, the 

states which influence the shape of a fusion barrier distribution. The analysis of barrier 

distributions thus offers a criterion for the relevance of the 'counter term' prescription in 

the Caldeira-Leggett approach. 

PACS numbers: 25.70.Jj, 24.10.Eq, 21.60.Ev, 27.60.+n 



Quantum tunneling plays an important role in a range of diverse phenomena in physics 
and chemistry. Recent attention has been focused on tunneling in systems with many 
degrees of freedom[||, ||. One of the interesting aspects of the problem is in determining 
which of the multitude of degrees of freedom must be explicitly included in any theoretical 
description, and which can be omitted. In particular, it is essential to define the role of 
excitation energy, or the degree of adiabaticity, in limiting the effectiveness of a specific 
degree of freedom. 

In nuclear physics, heavy-ion fusion reactions at energies near and below the Coulomb 
barrier provide an ideal opportunity to address this question. In order that fusion reac- 
tions occur, the Coulomb barrier created by the strong cancellation between the repulsive 
Coulomb force and the attractive nuclear interaction has to be overcome. Extensive exper- 
imental as well as theoretical studies have revealed that couplings of the relative motion 
of the colliding nuclei to nuclear intrinsic degrees of freedom cause enhancements of the 
fusion cross section at subbarrier energies, sometimes by several orders of magnitude, over 
the prediction of a single-barrier penetration model0. 

In a simple eigenchannel approach, such couplings result in the single fusion barrier 
being replaced by a distribution of potential barriers. A method of extracting barrier 
distributions directly from fusion excitation functions was proposedjQ, and stimulated 
precise measurements of the fusion cross sections for several systems. In a rigorous theo- 
retical interpretation, the barrier distribution representation is valid only if the excitation 
energy of the intrinsic motion is zero]^. Nonetheless, this method has been successfully 
applied to analyse data from heavy-ion collisions, where each intrinsic degree of freedom 
carries a finite excitation energy. These analyses of the barrier distributions have beauti- 
fully demonstrated the effects of coupling of the relative motion to surface vibrations[^, 0|, 
their multi-phonon excitations^, 0], deformations and the associated rotational excita- 



tions [0, [T0[ of the target nucleus and transfer reactions between the colhding nuclei||, 0]. 
Despite these successes, there are apparent conflicts regarding the role of projectile 
excitation. Each barrier distribution for the reactions ^''Ca + ^^^Pt, ^^^Os shows a charac- 
teristic structure, with a higher energy peak which has been associated with the octupole 
excitation of ^°Ca||ri|. Calculations of fusion cross-sections for the reactions ^^O + ^^^Sm, 



"^Ge in []12| , p!3| indicated that excitation of ^^O is important. In marked contrast, there are 
no specific features in the measured barrier distribution for the ^^O + ^^^Sm reaction]^ 
which can be associated with the excitation of ^^O; indeed it was shown in Ref. that a 
good theoretical representation of the barrier distribution is destroyed when the projectile 
excitation is included. 

All the above conclusions were based on comparison of the experimental results with 
simplified coupled-channels calculations. The simplification has been achieved by using 
one or more of the following approximations : 



1. the no-Coriolis approximation |1^ , where the centrifugal potential is assumed to be 



the same for all channels and equal to that in the elastic channel; 

2. the linear coupling approximation, where the nuclear coupling potential is assumed 
to be linear w.r.t. the coordinate of the nuclear vibrational excitation; 

3. the constant coupling approximation, where the coupling potential is assumed to be 
constant over the interaction range; 

4. intrinsic excitation energies are assumed to be negligible or are treated approxi- 
mately. 

The first approximation, common to most coupled channel calculations, including 



those presented in this Letter, has been shown to work well for heavy ion fusion calculations[15] 



Simplified coupled channel calculations]^, ^ |T3|, |T6| use the second approximation in con 



juction with either the third or fourth. Recent studies |T^, |18| have shown the linear 
coupling approximation is not valid even in systems where the coupling is weak, and 
that higher order couplings strongly influence the calculated barrier distributions. It is 
therefore probable that in reactions with nuclei like ^^O and ^°Ca, where the couplings 
to the octupole vibrational excitations are strong, barrier distributions calculated with 



simplified coupled-channel codes like CCFUS[|T^ do not provide a good representation of 
the fusion process. 

In this Letter we present the results of realistic coupled-channels calculations which 
demonstrate the effects of non-linear coupling and finite excitation energy of intrinsic 
nuclear (environmental) degrees of freedom, and resolve the apparently conflicting con- 
clusions regarding the influence of the projectile excitation. The relevance of the 'counter 
term' prescription of Caldeira and Leggett|l[| in heavy- ion fusion reactions is also discussed 
and the double counting problem of coupling effects is clarified. 

The coupled channels equations are solved by imposing the incoming wave boundary 
condition to simulate the strong absorption inside the fusion barrier. The real nuclear 
potential is assumed to have a Woods-Saxon shape and the depth was chosen to reproduce 
the experimental fusion cross sections at high energies using the single-barrier penetration 
model. The values of deformation parameters are extracted from the reduced transition 



probabilities using Px = 3^j;a y \,2 ■ The parameters of the calculations are listed in 
Table 1. 

In order to show the inadequacies of the often used linear coupling approximation, 
calculations were performed for the ^^O + ^^^Sm reaction using the linear coupling ap- 
proximation. The results of our calculations for the fusion excitation function and the 
barrier distribution are shown in Fig. 1. In the following discussion we concentrate on 



the latter since they are a more sensitive way to compare experimental data and calcu- 
lations. The dotted line shows the result when the excitation of ^^O is not included in 
the calculations. This calculation well reproduces the features of the experimental barrier 
distribution. Calculations including the excitation of the lowest-lying octuple state of ^^O 
are shown by the solid line. Even though the experimental barrier distribution around the 
lower energy peak (~ 60 MeV) is reproduced, significant strength is missing around the 
higher energy peak near 65 MeV. A similar discrepancy between theory and experimental 
data was encountered in Ref. 0, where calculations, shown by the long-dashed line, were 
performed using a modified version of the CCFUS code. It is clear that both calculations 
which include the octupole excitation of ^^O in the linear coupling approximation fail to 
reproduce the experimental barrier distribution. 

Realistic coupled channels calculations were then performed, where the couplings to 
the octupole vibrations of both ^^O and ^^''Sm are treated to all orders; i.e. the nuclear 
interaction is not expanded with respect to the deformation parameter, and the cou- 



pling matrix elements are numerically evaluated at each internuclear separation [|T^. It is 
remarkable that these calculations, shown in Fig. 2, re-establish the double-peaked struc- 
ture seen in the experimental data, which was absent in the equivalent linear coupling 
calculations. The shape of the barrier distribution obtained by including the octupole 
vibration of ^^O using all order coupling is now very similar to that obtained by ignoring 
it. This similarity becomes particularly evident when the calculated barrier distribution 
is shifted in energy, as shown by the dashed line in the figure. This is consistent with 
the general conclusion that the main effect of the coupling to inelastic channels whose 
excitation energies are larger than the curvature of the bare fusion barrier, i.e. an adia- 
batic coupling, is to introduce a static potential shift [|T^, ^, and hence, the shape of the 
barrier distribution does not change unless the coupling form factor itself has a strong 



radial dependence. 

In macroscopic quantum tunneling in condensed matter physics, the so-called counter 
term is often introduced in order to compensate for the static potential renormalization 
due to the coupling to the environment [0] . In contrast, in heavy-ion reactions, one usu- 
ally estimates the bare potential, for example by fitting the fusion cross section at high 
energies, and discusses the effects of channel coupling without introducing the counter 
term. Fig. 2 shows that this approach reproduces the experimental fusion cross sections 
and fusion barrier distributions without explicitly taking into account the excitation of 
the octupole vibrational state of ^^O. This indicates that the effects of its excitation is 
already included in the bare potential. If this is the case, the effect of the coupling to the 
3^ state of ^^O is double counted if the coupled channels calculations explicitly take it 
into account, resulting in a dramatic overestimate of the the experimental cross-sections. 
A recipe to cure this problem is to introduce the counter term as in condensed matter 
physics. Since the experimental data are well reproduced when the calculated distribu- 
tions are shifted to higher energies by 2 MeV, this shift evidently mimics the effects of 
the counter term. 

In contrast to the ^^0+^'^^Sm case, the analyses of ^°Ca + ^^^Os and ^^^Pt reac- 
tions, also performed using simplified coupled channel calculations []TT|, suggest that the 



excitation of ^°Ca is important in determining the observed barrier distribution. An im- 
portant difference between the ^^O and ^°Ca projectiles is that the excitation energy of 
the octupole vibration in the latter is smaller and nearly equal to the energy scale of the 
curvature of the fusion barrier, hence the coupling is intermediate between adiabatic and 
sudden. It is therefore interesting to investigate the degree of adiabaticity of the octupole 
excitation of the ^°Ca projectile. 

The results of the coupled channels calculations are compared with the experimental 



data in Fig. 3. All order couplings to both the target and the projectile excitations 
have been included. Although ^^^Pt and ^^^Os are transitional nuclei which lie between 
the 7-unstable and rotational limits in the interacting boson model[^, we have assumed 



that they are rigid rotors with axial symmetry. The ground state rotational band of the 
target nucleus, with states up to the 10^ member, has been included in the calculations. 
When the ^''Ca excitation is ignored, barrier distributions are obtained which are similar 
to those expected for a classically deformed nucleus and these are inconsistent with the 
experimental data. When the octupole excitation of ^°Ca is included in the calculations, 
a higher energy peak is introduced which agrees well with that observed in each reaction. 
The mutual excitation channels up to 4'^{^3~, the former and the latter refering to the 
targets and the projectile respectively, are also included in the calculations. It is apparent 
that the projectile excitation significantly affects the shape of the barrier distribution in 



this case, as suggested in the simplified coupled channel calculations in Ref. |[TT| . 

As has been shown in the discussions for ^^O + ^^^Sm reactions, the correct treatment 
of the coupling, without making the linear coupling approximation, significantly reduces 
the effect of projectile excitation on the shape of the barrier distribution. Calculations 
of the CCFUS-type, which fail in these regards, would therefore be expected to predict 
larger coupling effects than observed experimentally. The apparent success of the CC- 



FUS calculations reported in Ref. [|Tl|] was probably due to the compensation for this 
overestimate by the use of a smaller deformation parameter than that obtained from the 
octupole transition strength. 

The theoretical calculations for the reactions with the ^°Ca projectile still significantly 
underestimate the fusion cross section at low energies, even after the excitation of the 



projectile is taken into account. As suggested in Ref. [irT|, coupling to transfer channels 



which have been ignored in the present calculations, might enhance the fusion cross section 



at low energies. 

In summary, we have performed coupled-channels calculations for the fusion reactions 
^^O + ^^^Sm and "^°Ca + ^^^Pt, ^^^Os. The calculations with full order coupling show that 
the dominant effect of the excitation of the ^^O octupole state at 6.1 MeV is to renormalise 
the static potential without significantly changing the shape of the barrier distribution. 
On the other hand, the excitation of the 3" state at 3.7 MeV in '^'^Ca introduces well 
defined peaks in the barrier distribution. These results suggest a natural limit to the 
energy of states which need to be considered explicitly in coupled-channels calculations. 
The myriad of weak, high energy excitations which might be possible, contribute only to 
a potential renormalisation without affecting the shape of the barrier distribution. The 
effects of these excitations can then be included in the bare potential in coupled channels 
calculations. If these channels are explicitly included in the coupled channels calculations 
without introducing the counter term, they could be double counted depending on the 
choice of the bare potential. 

It has been shown that in order to interpret the high precision fusion excitation func- 
tions, that have recently become available, it is vital to perform exact coupled channel 
calculations which treat the excitation energy and the radial dependence of the coupling 
form factor correctly. Whilst CCFUS-based calculations have apparently been very suc- 
cessful in reproducing observed barrier distributions, it is clear from our results that care 
must be taken in their interpretation; the approximations used are unreliable even for rel- 
atively weak coupling strengths. Exact coupled-channels calculation is the only reliable 
means of quantitatively understanding the barrier distributions. 
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Table 1: Parameters used in the coupled channels calculations for the indicated reactions. 







Channel Couplings 


Potential parameters 


Reaction 


Nucleus 


Type A" E* (MeV) Px 


V (MeV) ro (fm) a (fm) 


1^0 + i^^Sm 


i^^Sm 


vib 3" 1.81 0.205 


105.1 1.1 0.75 




160 


vib 3" 6.13 0.733 




40Ca + i94pt 


40Ca 


vib 3- 3.70 0.339 


330.0 1.0 0.84 




194p^ 


rot 2+ 0.328 p2 = -0.154 

/34 = -0.045 




40Ca + 1920s 


1920s 


rot 2+ 0.206 P2 =0.167 

f]^ = -0.043 


148.0 1.1 0.84 



Figure Captions 

Fig.l: Fusion excitation functions (upper panel) and the barrier distributions (lower 
panel) for the ^^O + ^^^Sm reaction. The experimental data (filled circles) are taken from 
Ref. [0. The linear coupling approximation is used in the coupled channels calculations. 
In all calculations, the effects of the octupole vibration of ^^'^Sm are taken into account. 
The dotted line is the results when ^^O is treated as inert. The solid line is the result of 
the coupled channels calculations when the coupling to the octupole vibration of ^^O is 
also taken into account; the dashed line is the result of an equivalent CCFUS calculation. 
Fig. 2: As Fig.l, but for the case when the coupled-channels calculations have been 
performed by including all order coupling. The meaning of the solid and the dotted lines 
is the same as in Fig. 1, while the dashed line is the same calculation as the solid line 
with the average barrier increased by 2 MeV. 

Fig. 3: The comparison of the experimental fusion cross sections (upper panels) and fusion 
barrier distributions (lower panels) for the ^°Ca + ^^^p^^^ i920s reactions with the coupled 
channels calculations. In all calculations, the effects of the excitation of the target nuclei 
are treated in the rotational model and couplings to all orders are included. The dotted 
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lines are the results when '^''Ca is treated as inert. The solid lines include the coupling to 
the octupole vibrational state in ^°Ca. 
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